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Various possibilities are currently under discussion to explain the observed weakness of the intrinsic magnetic field of planet Mercury. One 
of the possible dynamo scenarios is a dynamo with feedback from the magnetosphere. Due to its weak magnetic field Mercury exhibits 
a small magnetosphere whose subsolar magnetopause distance is only about 1.7 Hermean radii. We consider the magnetic field due to 
magnetopause currents in the dynamo region. Since the external field of magnetospheric origin is antiparallel to the dipole component 
of the dynamo field, a negative feedback results. For an af2-dynamo two stationary solutions of such a feedback dynamo emerge, one 
with a weak and the other with a strong magnetic field. The question, however, is how these solutions can be realized. To address this 
problem, we discuss various scenarios for a simple dynamo model and the conditions under which a steady weak magnetic field can be 
reached. We find that the feedback mechanism quenches the overall field to a low value of about 100 to 150 nT if the dynamo is not 
driven too strongly. 
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1 Introduction 

The recent flybys of the MESSENGER spacecraft at planet Mercury confirm the existence of a large scale 



magnetic field (Anderson et al. 2009). The dipole surface field, however, is roughly one to two orders of 



theory ( 


Wicht et al. 


2007 Olson and 


igreement ( 


Heimpel et aZ.|2005 


Stanley 


Glassmeier et al. 


200 


7) with different 



dynamo configurations. Here, we further study the feedback dynamo scenario suggested by Glassmeier 



et al. (2007) who investigated the interaction of the dynamo and the magnetospheric field. They derived 
two stationary solutions and ascribed the weaker solution to Mercury's magnetic field. They however do not 
address the question how the dynamo reaches either of these solutions. Allowing a variable magnetopause 
which depends on the internal field and solar wind conditions, it is so far not conceivable how a dynamo 
can develope into a state where it can be quenched by the external feedback field. Therefore, the present 
study aims at discussing conditions under which a steady and weak magnetic field can evolve when the 
dynamo is exposed to a magnetospheric magnetic field. 
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Magnetopause 

Figure 1. Schematic illustration of the feedback mechanism. The planet's dynamo generates an internal field £?j n t. The interaction with 
the solar wind causes a magnetopa use current Jmp which itself induces an external field B ex t which is of opposite orientation to the 
internal field in the Hermean core jGlassmeier et al. |2007[ |. 

2 A Hermean feedback dynamo 

The magnetopause currents caused by the interaction of Mercury's magnetic field with the solar wind 
generate an external field which reaches into the planet's interior. Since the internal magnetic field is weak, 
the magnetopause is located close to the planet. We thus expect a stronger external field contribution 
in the dynamo region than, for example, for Earth. In the terrestrial case, the subsolar magnetopause 
is located at about 10 planetary radii and its influence on the internal dynamo process is negligible. In 
contrast, at Mercury the magnetopause is close to the planet and the external field has to be taken into 
account in the solution of the dynamo problem. 

The relative orientation of the internal and external magnetic fields is of significance. As seen in figure 
[TJ the magnetopause currents generate a field canceling the field outside the magnetosphere. Inside the 
magnetopause internal and external fields are parallel. At the core-mantle boundary the situation is 
different. The internal dipole field possesses a vector component along the rotational axis of the planet 
that is anti-parallel to the externally generated magnetic field. Thus, a negative feedback situation results. 

Since the feedback field is stronger for a close magnetopause, we concentrate on a weak initial dynamo 
field. This situation corresponds to the onset of dynamo action or to the time period after a polarity 
reversal when the dipole field is weak compared to higher multipoles. In order to gain first insights into 
the system's possible temporal evolution, we reduce its complexity by coupling a simple kinematic internal 
dynamo to an idealized external magnetospheric field. 



3 Response function 



The external field arising from magnetopause currents depends on the distance of the magnetopause to the 
planet and the spatial current distribution. The dynamical magnetopause position, parameterized by the 
stand-off distance R s at the subsolar point, is mainly determined by the pressure equilibrium between the 
planetary magnetic pressure, with the dipolar part as the main contribution and the solar wind dynamic 



pressure (e.g. Baumjohann and Treumann 1996): 



Rm 



'2(9?,ir 



1/6 



MPs 



(1) 



April 30, 2010 1:30 Geophysical and Astrophysical Fluid Dynamics GAF'alpha'omega'v5r4 



The initial temporal evolution of a feedback dynamo for Mercury 3 

Here Rm, Mo, 9\. int an d Psw denote the Hermean planetary radius, the permeability of free space, 
the internal axial dipole Gauss coefficient and the solar wind ram pressure, respectively. Equation ([I]) 
demonstrates that the stand-off distance depends on the internal field strength like (^i int) 1 ^ 3 - -^ ne 
magnetopause is thus located close to the planet for weak magnetic fields like the one found at Mercury. 
In contrast to that, Earth exhibits a relatively strong magnetic field with a distant magnetopause and 
negligible influence on the internal dynamics in the planet's core. 

In general, when the shape of the magnetopause and the planetary dipole field strength are given and 
the solar-wind is assumed to be field-free, the external field from magnetopause currents can be calculated 
without explicitly determining the currents. This is achieved by shielding the internal field by an external 
potential field at the magnetopause, such that the magnetic flux through the magnetopause vanishes. 
The field-free approximation is applied since incorporating the ever-changing interplanetary magnetic 
field (IMF) characteristics would require detailed hybrid modeling of the solar wind interacting with the 
planetary magnetic field or long-term in-situ magnetic field observations which are not available at this 
time. Altogether, in order to determine the external field the stand-off distance must be known and then 
the internal dipole field strength sets the shielding current strength. It is therefore possible to express the 
external field strength B cy ± as a response function / of the internal one maintained by the dynamo process: 



B 



/«int) 



(2) 



There exist several models for the external field for various solar wind and planetary magnetic field 



conditions. For the Hermean case, Glassmeier et al. (2007) constructed a simple model treating the 



magnetopause current as a circular line current in the equatorial plane. The well-studied terrestrial 



situation can be described with a semi-empirical model by Tsyganenko and Sitnov (2005). In that study, 



the stand-off distance and the current strength depend on solar wind conditions. Making use of the 
aforementioned field-free approximation, the model prescribes the magnetopause shape as an ellipsoid 
with a cylindrical continuation as a magnetotail. The spatial parameters of this magnetopause are fitted 
to satellite observations. At this boundary the magnetic field of the planet represented by its dipolar part 
and contributions arising from several magnetospheric current systems are partially shielded depending 
on IMF conditions. 



For a more realistic representation than Glassmeier et al. (2007) we adapt the terrestrial model of Tsy- 



ganenko to Hermean conditions following the approach of Korth et al. (2004). First, we assume a centered 



axial dipole as the planet's intrinsic magnetic field. Furthermore, since there are no permanently trapped 
particles expected because of the low internal field strength of Mercury, we neglect the magnetospheric 
ring current. In order to scale this Tsyganenko model to Hermean conditions we make use of the scaling 
law 



B M (r M ) = B E (r E ) = B E (Br M ) with k = k p k b 



Psw,m\°- 14 ( 30,574nT 



1/3 



2nPa/ 



yi,int,M 



(3) 



Here Bm^m), Be(te) and g{ ktM denote a magnetic field vector in the Hermean system at the position 
tm, a magnetic field vector in the terrestrial system at the position te and the internal, axisymmetric 
dipole Gauss-coefficient of the magnetic potential expansion. The scaling factor k p due to different solar 



wind ram pressures at Earth and Mercury has been extrapolated from a model by Tsyganenko (1996) 



fitting satellite observations at different solar wind conditions to their model magnetopause. At Mercury 
the solar wind ram pressure is taken to be 13.4 nPa assuming an average solar wind speed of 400 km/s 
and an average proton number density of 5 x 10 7 m -3 (e.g. Glassmeier 1997). The value of 2 nPa is the 



average solar wind ram pressure at Earth. The second factor Kb, accounting for the different magnetic 



moments of the planets, has already been used by Korth et al. (2004). The factor 30,574 nT is the 



terrestrial magnetic dipole field strength at the equator around 1980 as it has been used in the Tsyganenko 



(1996) model. The linear factor kb thus scales the magnetospheres in such a way that the equatorial field 
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Figure 2. External dipole field strength as a function of the internal dipole field strength. The crosses indicate the values computed with 
the scaled Tsyganenko model and the solid line marks the model function with incorporated transition. The greyed area visualizes this 
transition region where the mode of interaction between the solar wind and planet is changed due to the planetary surface. 



strengths in both planetary systems are equal. 



The described magnetospheric model provides the full spectral multipole information of the external field 
in response to any internal dipole field strength. The dynamo, however, is only affected by a long-term 
average magnetopause field. The time span needed for an external field to diffuse through the entire core 
region is of the order of 



I? 
V 



fi aL 2 Ri 35 x 10 3 a , 



(4) 



where r\ is the magnetic diffusivity, cr = 6xl0 5 S/mis the assumed electrical conductivity (Suess and 
Goldstein 1979 ) and L = R cr ab — Rich is the radius of the outer core shell with i? cm b = 1860 km as the 
radius of the core-mantle boundary ( Spohn et aZ.||200T ) . The inner core radius i?; c b is not well constrained 
but an Earth-like value of 0.35i? cm b is chosen here. Therefore, any external non-axisymmetric magnetic 
multipole contribution to the overall field in the planet's interior would cancel over the Hermean orbital 
rotation period (88 Earth-days) and planetary rotation period (59 Earth-days). 



As external multipoles of degree I > 1 decrease towards the planet, we furthermore restrict ourselves for 
simplicity reasons to the strongest multipole 1 = 1. This provides a uniform magnetic field in the interior 
which is aligned with the planetary rotation axis. As an estimate of this external field we take the magnetic 
field value obtained from the model at the subsolar point on the Hermean equator. The resulting hyperbolic 
response function is shown by crosses in figure [2j where the external and internal fields are represented by 
their multipole coefficients g^ ext and g\ 

inf respectively. For a strong internal field, the magnetopause is 
pushed further away from the planet, thus resulting in a weak external feedback. In contrast, the external 
feedback is strong for weak internal fields since the magnetopause is closer to the planetary surface. 
We modify the response function to exclude the unrealistic case of a stand-off distance located within 
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the planet, R s < R M , which is equivalent to g\ int < 91.8 nT according to 



take into account the finite extent of the magnetopause (Berchem and Russel 



1|). Furthermore, we need to 
982 ) . This implies that the 



magnetopause currents are distributed over a finite radial extent. Some of the consequences for magnetic 
field measurements in the Hermean system are discussed by Glassmeier et al. (2010). While detailed 



modeling of the complex magnetopause current structure is beyond the scope of the present study, we take 
a first step to respect the finite thickness. The full response function is assumed to apply only at distances 
greater than 500 km away from the planetary surface. This corresponds to a typical magnetopause thickness 
( Berchem and Russell|1982 ). In consequence, the planetary dipole coefficient must exceed about 161 nT to 
yield an entirely undisturbed magnetopause. For weaker fields, i.e., for smaller R s , we modify the response 
function such that a smooth transition towards the R s < Ryi situation emerges. The solid line in figure 
[2] shows this modified feedback function, whose functional form is described in the following. Throughout 
the transition interval 91.8 nT < fl? int < 161 nT we model the external field with a response function 



5i°ext = 9-00 x 10 6 (g° 1M - 91.8 nT) (<7?, int /nT) 



-2.73 



(5) 



while for the remaining interval g® illt > 161 nT the feedback function is parameterized with an exponential 
function 



yi,ext 



1.37 x 10 5 nT exp (-5.96 (g? int /nT) 



^5.27x10- 



(6) 



fitted to the findings from the Tsyganenko model. This parameterized response function allows us to 
calculate the influence of the magnetospheric magnetic field on the dynamo without explicitly evaluating 
the stand-off distance and the modified Tsyganenko model. 



4 An aJl-dynamo embedded in an external magnetic field 



In order to describe the influence of an imposed external magnetic field on the dynamo process, an addi- 



tional induction term in the dynamo equation is introduced (Levy 1979, Glassmeier et al. 2007): 



dB 

~dt 



V x [v x (B + B ext )] + ?? AB 



(7) 



where v denotes the velocity, 77 the magnetic diffusivity, B the dynamo field and B ext the external 
magnetic field. To study the temporal evolution of the feedback dynamo we adapt a version of a 
ID kinematic mean-field a£l- model presented by Schmitt and Schussler (1989), who studied different 



non-linear quenching mechanisms with application to the Sun. With the magnetospheric feedback we 
introduce another non- linear quenching method but within the context of a planet with a magnetosphere. 
The main scope of this paper is to address the question how the coupled system can dynamically evolve 
into a weak field solution. The simple kinematic dynamo serves to get a first picture of the various 
scenarios that may arise. 



The model considers dynamo action in a differentially rotating spherical shell with an outer core radius 
of R C mb- The radial variation of the magnetic field and of the induction effects are neglected. About the 
latter little is known in the case of Mercury and any specification seems arbitrary. The neglect of a radial 
dependence of the magnetic field is only permitted for a thick shell, as it is probably the case for Mercury's 
fluid core. Furthermore, we assume rotational symmetry, so that all quantities are independent of the 
azimuth (d v = 0) and thus depend solely on the colatitude 9. The magnetic field is decomposed into a 
poloidal component, described by a vector potential A = (0, 0, A) and a toroidal magnetic field component 
(0,0, B) by 



B = {0,0, B) + V x (0,0, A) 



(8) 
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The toroidal field is produced by a constant radial shear d r Vt = Of = const, through the so-called Sl-effect. 
The poloidal field is maintained by the a-effect, a parameterized interaction of small-scale field and 
small-scale flow. We assume a simple harmonic dependence a(9) = oqcos^. 



In order to compute the magnetic field in the magnetosphere, the field must be continued outside the 
dynamo shell. Since the radial component of the magnetic field must be continuous at the core-mantle 
boundary we can analyze B r = (V x A) r to find the internal dipole Gauss coefficient <7i int . Any possible 
influence of the embedding electrically conducting mantle is neglected here. With the internal dipole 
coefficient known at each time step we computed the magnetospheric response with the parameterized 
response function for the successive time step. The equations are non-dimensionalized by using the magnetic 
diffusion time scale r = R 2 m b/v^ the length-scale -R cm b and an appropriately chosen magnetic scale Bq. 
Furthermore, we abbreviate A{&) = A(9)sm.9 and B(9) = B(9)s'm9. For the dimensionless uniform 
external field we choose 



iext(0) 



ffl,ext 

~Bo" 



COS O 8111 ( 



0) 



The ail-dynamo equations with an ambient poloidal magnetic field following the approach by Levy ( 1979 ) 
are written as 



d t A = d 2 6 A - cot 9 d e A + cos 9B 
d t B = dlB - cot 9 d e B + Psinf? ( d e A - A e 



(10) 
(11) 

with the poloidal external field contributing to the induction effect acting on the toroidal component. 



The first two terms of the right-hand side of equations (10) and (11) describe the diffusion of the 



poloidal and toroidal field, respectively, the third term of (10) the action of the a-effect on the toroidal 



field and the third term of ( 11 ) the differential rotation acting on the internal and the external poloidal field. 



The model is controlled by the dimensionless dynamo number P = R^^Q' oq / rj 2 . Without an external 
field the magnetic field grows when P exceeds a critical value of P CT it = 46 and decays otherwise. At 
values P S> -Pcrit the dynamo would begin to show an oscillating behavior. The complete mode structure 



is described by Schmitt and Schiissler (1989). It is qualitatively also typical for a 2D thick layer dynamo 



see e.g. Parker 1971). In the present study we are interested in the monotonically evolving mode for 



dynamo numbers P cr i t < P < 70. 

A toroidal magnetic field of the form -B see d = 10~ 5 sin 9 with a low amplitude compared to Bq serves 
as seed field in order to avoid a dependency of the results on the properties of the initial field. Temporal 



integration of the modified induction equations (10) and (11) is computed numerically using a finite 



differencing scheme, where the diffusion terms are treated implicitly and the induction terms explicitly. 
The temporal evolution for different P with and without feedback is presented in figure [3j In a fully self- 
consistent dynamo the growth of the magnetic field is limited by the Lorentz force acting back on the flow. 
Since this feedback mechanism is missing in this kinematic model, supercritical dynamo numbers P > P C nt 
result in an unbounded exponentially increasing magnetic field strength over time as exemplified for the 



two cases in figure |3j For example, Schmitt and Schiissler ( 1989 ) limit this growth by an a-quenching 
whose effect is shown as the dashed line in figure |3| 

The negative feedback from the external field provides an alternative quenching mechanism. We generally 
start with a small internal field which is insufficient to produce a magnetopause above the surface and thus 
provides no quenching. For supercritical dynamo numbers the internal field then grows until the external 
field has developed sufficiently to provide the necessary quenching as visualized in figure [3j This results 
in a stationary solution with a magnetic field strength that depends on the dynamo number. This can 
be seen comparing the P = 50 and P = 54.5 cases in figure |4j The first one is stabilized after about 5 
diffusion times, the latter after about 25 diffusion times. The saturation level is between 100 and 145 nT. 
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Figure 3. Comparison of the temporal evolution of the internal dipole strength depending on different dynamo numbers P and feedback 
turned on or off. 
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Figure 4. Temporal evolution of the internal dipole strength for various dynamo numbers P with the magnetospheric feedback turned 
on. 



However, when the dynamo number is chosen bigger than P = 54.5 the quenching is insufficient and the 
exponential growth is only delayed. This happens when the internal dipole strength exceeds <?i mt = 145 
nT where the external field reaches its maximum. This level is marked in figure [4] with a dotted, horizontal 
line. We note that the starting field strength has to be lower than 145 nT for the quenching to work. The 
saturated field strength is independent of the initial amplitude. The duration of the delay for P > 54.5 and 
the ultimate exponential growth rate both depend on the dynamo number. At P = 54.7 the exponential 
growth phase sets in after about 60 diffusion times whereas the delay has virtually vanished at a dynamo 
number of P > 70. 

The range of dynamo numbers for which the magnetospheric field limits the growth of the dynamo, 
here the relatively small parameter range of 46 < P < 54.5, depends strongly on the maximum value 



April 30, 2010 1:30 Geophysical and Astrophysical Fluid Dynamics GAF'alpha'omega'v5r4 



Heyner et al. 



120 
100 

80 

60 

40 

20 





1000 



100 * 



100 200 300 400 500 
9?. im [ nT ] 



70 



69.5 



-69- 



P=60 



20 40 60 80 100 
Time [diffusion time] 





140 


■ i i | i i i | i i i | i i i | i i i 

51 




130 




I — 

c 


1 PO 

1 C\J 




c 




/ / model 3: P=50 


XJi 


110 


I 51 




100 






l model 1 : P=50 




90 





20 40 60 80 100 
Time [diffusion time] 



Figure 5. Parameter study with artificially altered response functions. Left panel: Model 1 corresponds to the already used form for 
Mercury. In model 2 the function of model 1 is amplified by a factor of 2. In model 3 the function has the same amplitude as in model 
1 but the maximum is shifted to a higher internal field value. Middle panel: Temporal evolution of the internal dipole field strength for 
various dynamo numbers P for model 2. Compare with figure 4 for model 1. Right panel: Stationary saturation levels of the internal 
dipole strength for two dynamo numbers P for model 3 in comparison with model 1. 



of the response function. In order to analyze the dependence of the response function, we investigate 
two different cases with modified response functions. If we multiply the response function by a fac- 
tor of 2 (model 2 in figure [5]), dynamo numbers up to 69 result in stationary solutions as is shown in figure [5] 

The field strength at which stationary solutions saturate, in the original model (hereafter referred to as 
model 1), values between 100 and 145 nT, is determined by the rising part of the response function. If 
the maximum is at higher values of g± int (model 3 in figure pi), a higher saturation field strength results 
compared to model 1 for the same dynamo numbers. 



5 Conclusion and outlook 



Using a kinematic afi-dynamo in a feedback configuration, we have demonstrated that the feedback of 
the external field on the internal dynamo mechanism can indeed result in relatively small field strengths 



below 150 nT as suggested by Glassmeier et al. (2007). However, in our simplified kinematic dynamo 



model the responsible quenching would only be sufficient in a narrow regime where the dynamo number 
does not exceed 18% of its critical value. If Mercury is captured in the quenched regime our model 
implies that the Hermean dynamo is unique. It should be noted here that alternative explanations for the 



weak Hermean dynamo field (e.g. Stanley et al. 2005, Christensen 2006) also require the assumption of 
special conditions for Mercury. The saturation field strength strongly depends on the assumed response 
function describing the dependence of the external field on the internal field strength. Unfortunately, 
very little is known about the underlying interaction, especially for a magnetopause close to the surface 
which would be appropriate for Mercury which is neccessary for our suggested feedback mechanism to work. 



This paper is part of a series of studies examining the model of a feedback dynamo scenario. Glassmeier 



et al. (2007) made use of extensively simplified models and examined stationary dynamo solutions without 



addressing the question how these stationary solutions could be realized. This problem has been addressed 
in this study. We further consider an analytical solution to an approximation of the kinematic dynamo 
problem which allows us to examine the influence of the shape of the response function on the dynamo 
solution. The results could be useful for the application of the idea of a feedback dynamo to other astro- 
physical bodies such as gas giants close to their host star. Furthermore, we address the response function 
(also for higher magnetic multipoles) for Mercury by using a hybrid code simulating the interaction of 
Mercury's magnetosphere with the solar wind. Another investigation concerns how a three-dimensional, 



self-consistent, numerical dynamo model in approximate magnetostrophic balance (Wicht 2002) reacts 
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to an imposed uniform and constant-in-time external field. From the results of these simulations we will 
know what kind of characteristic reactions of the dynamo we can expect when examining the full time 
dependent, 3D model with the exact and full magnetospheric response function. 
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